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Abstract 

Summation by parts is used to find the sum of a finite series of generalized 
harmonic numbers involving a specific polynomial or rational function. The 
Euler-Maclaurin formula for sums of powers is used to find the sums of some 
finite series of generalized harmonic numbers involving nonnegative integer pow- 
ers, which can be used to evaluate the sums of the finite series of generalized 
harmonic numbers involving polynomials. Many examples and a computer pro- 
gram are provided. 

Keywords: generalized harmonic number. 
MSC 2010: 11B99 

1 Harmonic Numbers and Sums of Powers 

Let the generalized harmonic number for nonnegative n, complex order m and complex 
offset c be defined as: 

" 1 

H ( c ™ } = Y 7 rr- (1-1) 

For ease of notation the following convention is used: 

= Hjg (1.2) 
The traditional harmonic numbers are a special case: 

H n = tf« (1.3) 

The following sum of powers can be expressed as a generalized harmonic number for 
nonnegative integer p: 

n 

£V = (L4) 

fc=l 

Let the be the Bernoulli numbers with =1/2 instead of B\ = —1/2: 

B+ = (-l) n B n = B n + 6 nA (1.5) 



1 



Then application of the Euler-Maclaurin summation formula 8J to the sum (1.4 1 
yields for nonnegativc integer p: 



H (-p) 



, p+i 

p) = 1 
p + 1 ^ 

1 k=l 



k i^-fe+i" 



(1.6) 



This formula expresses any power sum ( 1.4 1 with nonnegative integer p in a polynomial 



in n of degree p+1, and is used below to obtain a similar reduction for some finite series 
of generalized harmonic numbers. For a list of the resulting expressions of formula 
(l.Gl up to p = 5, see section [5] below. In all following derivations it is assumed that 

0° = 1 EflE]. 



2 Summation by Parts 

For application in section [3] below summation by parts is defined and proved. 

Theorem 2.1. Let a < b be two integers and let {xk} and {yk} be sequences of 
complex numbers. Let {s^} be the sequence of complex numbers defined by: 



Sk 



(2.1) 



Then there is the following summation by parts formula: 

6-1 6-1 



^2 XkVk = Sb-iy b - ^2 s k (y k+1 - y k ) 



(2.2) 



Proof. From (2.1) follows 



Xk = Sk - s k -i 



Then (2.2) becomes 



6-1 



Sb 



-lVb = ^2i(sk - s fc _i)y k + s k (y k+1 - y k )} 



k—a 
6-1 



^2(s k y k+1 - s k -iVk) 

k—a 

Sb-lVb - Sa-lVa 



(2.3) 



(2.4) 



Since s a _i = by definition (2.1 ), the theorem is proved 



□ 
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3 Generalized Harmonic Number Identities 



The following identity is a finite series of generalized harmonic numbers involving a 
certain polynomial or rational function. 

Theorem 3.1. For nonnegative integer n and complex m, w: 

n 

Y^[{k + 1)™ - k w ]H { k m) = (n + 1) W H^> - H^' w) (3.1) 

k=0 



(k + l) w — k w and = if! , where the Sk become: 



Proof. The summation by parts theorem (|2.2|) is applied with a = 0, b = n + 1, Xj, 

k 

Sfc = ^[(i+ir-r] = (fc+ir (3.2) 

i=0 



and (2.2) becomes: 



+ 1)- - k w ]Ht ] = (n+ irH™ - J> + l) w {k + 1)- 

k=0 k=Q 

-{n+1) H n+1 - H n+1 
= {n+ l) w H^ - H^ 1 -^ 



(3.3) 



□ 



Let the following sum be given [31 HI [SI E] : 

Y,\H k = l -{Hl + HW) (3.4) 

fc=i 

Then applying this theorem for m = 1 and w = — 1 and = Hn+\ ~ V( n + 1)™ 
yields: 

tjk!i flt = l^-4) (3 - 5) 

fc=0 

The following identities are finite series of generalized harmonic numbers involving 
nonnegative integer powers, which can be used to evaluate finite series of generalized 
harmonic numbers involving polynomials [51 [7] . 

Theorem 3.2. For nonnegative integer n, p and complex m: 

n 

^k p H ( k m) =F{n,p,m) 

(3-6) 

i p+ / _i_ i\ 

- JT(-P) JjM + Jj(m-p) — VP IB+ iJ 



fc=l 



(m — k) 
p—k+l ±J -n 
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Proof. The order of summation over k and i is reversed 

k 



n n 



Yk?Hl m) = Yk>>Y — = Y — 

Z ^ ^ Z ^ Z ^ j^m Z ^ Z ^ 

fc— 1 fc— 1 i— 1 2—1 fc— z 

The sum over fc is splitted into three parts: 



Y k p = y kP + iP - E kP 



(3.7) 



(3.8) 



k=i k=l 



k=l 



The last sum is expressed by ( 1.6 1, and the order of summation over k and i is reversed 
again. □ 



For a list of the resulting expressions of formula (3.6) up to p — 5 and m = 4 see 
section [U below. 

Theorem 3.3. For nonnegative integer n, p and complex m: 

n 

Y k p H^ k = G(n,p, m) = H^H™ + S^H™ 



k=0 



V 



k=l 



(3.9) 



[£+_ fc+1 + (P k + l)Hi k -^}Ht~ k) 



Proof. The order of summation over k and % is reversed 

n n~k 



n ^ 7i— % 



fe=l 



fc=l i=l 



1=1 fc=l 



The sum over k is evaluated with (1.6): 



n—i p+1 

w = ^—Y 



k 



k=l L " ' ' k=l 

On the right side the power is evaluated with the binomial theorem: 



(n-i)* = £(-!)' 



1=0 



n k ~ l i l 



The trinomial revision of binomial coefficients is used: 

'p+l\fk\ fp+l\fp-l + l 



I 



k-l 



(3.10) 



(3.H) 



(3.12) 



(3.13) 



The terms with I — are splitted off, the order of summation over k and I is reversed, 
and k and / are interchanged: 



Yw = — Y 



1 P+1 / I 1\ P+1 / 7 I 1 \ 



k=l 



k=l 



l=k 



(3.14) 



The first term is evaluated with (1.6 1, and for the second term the summation range 



over I is changed and (1.6) is used 

p-fc+i 



The order of summation over k and i is reversed again. 



(3.15) 



□ 



For a list of the resulting expressions of formula (3.9) up to p — 5 and m — 3 see 
section [5] below. 

4 Generalized Harmonic Numbers Identities with 
Nonzero Offsets 

The formulas above can be used to derive identities for generalized harmonic numbers 



(1.1) with nonnegative integer offset c = s. For this the following equation is used 

(4.1) 



with nonnegative integer s: 



If(m) _ fj(m) _ Tr(rn) 
n s,n — n s+n n s 



which leads to: 



fc=0 



k=0 



(4.2) 



= -H^(H^ + 5piQ ) + £ k PH s+k 



(m) 



fc=0 



The last sum can be evaluated as follows: 

n+s 



k=0 



k—s 
n+s 

= £(*-*)*# 



(m) 
fe 



s-1 



(4.3) 



PH r(m) 



-5> -*)*/?<*"> 



fc=0 



fe=0 



The power can be expanded with the binomial theorem, and changing the order of 



summation results with F(n,p,m) from (3.6) in: 

£ kPH *+l = E(-!) fc U K + - fc ' m ) - F(* - 1,P - fc, m)] (4.4) 

k=0 n ^ ' 



fc=0 



which can be substituted in (4.2). In the case s — the only nonzero term in this 
formula is for k = with s fc = U = 1, and because F(—l,p, m) = the result in this 
case is F(n 7 p,m) as expected. 
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For a list of the resulting expressions of formula ( |4.4[ ) with s = n up to p — 5 and 
to = 2 and with s = 2n up to p = 5 and m = 1 see section [5] below. 

Similar to the derivation above the second type of sums can be evaluated: 



n 



k=0 



fc=0 



(4.5) 



k=0 



The last sum can be evaluated as follows: 

n n+s 



s-1 



E fcp <L* = e kPH( s :i- k B» + 1 + *) p ^-U 



(4.6) 



fc=0 fc=0 fc=0 

The power can be expanded with the binomial theorem, and changing the order of 
summation results with G(n,p,m) from (3.9| in: 



" p / \ 

E kPH i+l-k = G ( n + ™) - E ( I) (" + !)^"G( S - 1, fc, 777) 

fe=0 k=0 ^ ' 



(4.7) 



which can be substituted in (4.5 1. 

For a list of the resulting expressions of formula (4.7) with s = n up to p — 5 and 
to = 1 see section [5] below. 



5 Examples 

Using the computer program of section [6] below, a list of many resulting expressions 
of the main formulas is provided. 



Formula ( 1.6 1 



= 77 (5.1) 

H^ = ^77(77+1) (5.2) 

ij(- 2 > = -77(77 + 1) (277+ 1) (5.3) 
6 

H^ = \n*{n + lf (5.4) 

i?, ( r 4) = 7^77(77 + l)(2n + 1)(3t7 2 + 377 - 1) (5.5) 

= ^n 2 (n+ 1) 2 (2t7 2 + 2n - 1) (5.6) 
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Formula (3.6 1 



fe=0 



H k = (n + l)H n+1 - (n + 1) 

fc=0 

" 1 
5] fci/fe = Ht 1] H n+1 - -n(n + 1) 

" 1 
53 k 2 H k = H^H n+1 - — n(n + l)(4n + 5) 

fc=0 

" 1 
^ fc 3 ff fc = Ht 3) H n+1 - —n(n + l)(n + 2)(3n + 1) 

n 1 

^ fc 4 tf fe = H^ 4 >H n+1 - j^n{n + l)(72n 3 + 243n 2 + 167n - 32) 

fe=0 

n 1 
53 fc 5 JJ fc = H^H n+1 - — n(n + l)(n + 2)(2n + l)(10n 2 + 19n - 9) 

f^fff = ^ + 1)^-^+! 
fc=0 

53 fe 2 iff = i^" 2 ^ - ±2^1 - i(n + l)(n - 1) 

fc=0 

53 fc 3 fff = ffH^+i - -n(n + l)(2n + 1) 

fe=0 

" 1 1 

53 fc 4 < > = /4- 4) #fi + ^^n+i - e> + + 2)(3n2 - n + 1} 

) 

n , 
53 fc 5 ijf = Jf^/f^ - — n(n + l)(12n 3 + 33n 2 + 7n 7) 

±H^ = (n + l)Hi%-H^ 

k=0 

n 

53 kH^ 3) = H^H^I-l + 2^+1 _ 2 Hn+1 

k=0 

± = H^H^l - \h% + \ Hn+l - i(„ + 1) 



fc=0 



fc=0 



(5.7) 
(5.8) 
(5.9) 
(5.10) 
(5.11) 
(5.12) 
(5.13) 
(5.14) 
(5.15) 
(5.16) 
(5.17) 
(5.18) 
(5.19) 
(5.20) 
(5.21) 



fc=0 



fe=0 



n 1 

k 3 H^ = HtVH^ \H n+1 - ±(» + l)(n - 2) 

fc=0 

n 1 1 

x: fc 4 4 3) - ^- 4) ^i + - ^ + i)(4n2 - n + 2) 

n 

2 fc 5 4 3) = fli- 5) ^Si + Y2 Hn+1 24 ( " + + 2)(n2 _ " + 1} 

n 

$>f = (n + l)H^ 1 -H^ 1 

k=0 

n 1"1 
fc=0 

J.2 rr(4) _ rr(-2) rr(4) 1 rr(3) , 1 rr(2) 1 „ 

2_^k H k — 'H n+1 - g-^n+l + 2^"+! - 3 " +1 

fc=o 

£ = - i ^ + ±If n+1 - i(n + 1) 

fe=0 

£ fc 4 < } = FJM)^ + iflW - - i(n + l)(n - 3) 

£ fe 5 < = H^H^% + l^x - ^(n + l)(2n 2 - 2n + 3) 



fe=0 



fc=0 



Formula (3.9| 



^ if„_ fe = (n + l)H n+1 - (n + 1) 

fc=0 

n „ 

^ fci/„_ fe = H^H n+1 - -n(n + 1) 

fe=0 

n 

J] fc 2 ff„_ fe = Ht 2) H n+1 - ~n(n + l)(22n + 5) 

fc=0 

" 1 
^ A: 3 i/ n _ fc = fT(- 3 )fl n+1 - — n(n + l)(25n 2 + 13n - 2) 

n 1 

V fc 4 i? n _ fe = H^H n+1 - — n(n + l)(822n 3 + 693n 2 - 133n - 32) 

lolIU 

fc=0 

n 1 

2J fc 5 #„-fc = H { - b) H n+l - ^nfa + l)(98n 4 + 116n 3 - 21n 2 - 19n + 6) 



(5.22) 
(5.23) 
(5.24) 
(5.25) 
(5.26) 
(5.27) 
(5.28) 
(5.29) 
(5.30) 

(5.31) 
(5.32) 
(5.33) 
(5.34) 
(5.35) 
(5.36) 



fc=0 



E^ fc = (« + l)^Si--ffn + i (5.37) 

fc=0 

" 1 1 

fe=0 

E fc 2 ^ fc = Ht^H^ - -(6n 2 + 6n + l)H n+1 + -(n+ l)(5n + 1) (5.39) 
E * 3 ^ fc = ^ ( r 3) ^Si - 3^- 2) ^„ + i + -n(n + l)(2n + 1) (5.40) 



fe=0 



fe=0 



fe=o JU (5.41) 
+ ^(n + l)(77n 3 + 65n 2 + n - 2) 

n 

fc=o (5.42) 
+ ^™(™ + l)(522n 3 + 633n 2 + 37n - 67) 

oDU 

n 

E ^n-fc = (« + - ( 5 - 43 ) 

fc=0 

E*^fc = - ^(2n + l)^ 2 i + ±ff n+1 (5.44) 

fc=0 

^fc 2 ^ fe -^ 2) ^- 6 (6n 2 + 6n + l)^ 1 + -(2n+l)i?„ +1 --(n + ^ (5.45) 



fe=0 



n 1 

E fc3 ^ fc = **- 8) *2!i - 3^- 2) 41 + V a + 6n + l)H n+1 

fc=o * (5.46) 

-i(„+l)(7n + 2) 



1 

E fc 4 ^* = H^Hi% -W + 60n 3 + 30n 2 - 1)^ 

fc=o M (5.47) 

+ n(n + l)(2n + l)#„+i - ^-(n + l)(94n 2 + 59n + 2) 

oU 

n 

E fc5 ^ = ^- 5) ^i- 5 ^" 4) ^i 

fe=o (5.48) 
+ ^ (30n 4 + 60n 3 + 30n 2 - l)#„+i - ^ (n + 1) (57n 3 + 57n 2 + 5n - 2) 



Formula (4.4) 



H n+k = (2n + l)H 2n+1 - nH n - (n + 1) (5.49) 

fc=0 

n 1 
J2 kH n+k = H^H n + -n(n + 1) (5.50) 

fe=0 

n l 

k 2 H n+k = Hl-^[2H 2n+1 - H n ] - —n(n + l)(10n + 11) (5.51) 

fe=0 

n 1 

V fc 3 iT n+fe = H^H n + —n(n + l)(7n 2 + 7n-2) (5.52) 
. — ' 48 

fc=0 

™ 1 

Y^k 4 H n+k = i?(- 4 )[2i7 2li+1 - J ff n ]-— — n^ + l)(282n 3 + 603n 2 + 2577i-92) (5.53) 

fc=0 

n 1 

5^ fc 5 iJ n+fc = Ht 5) H n + ^n(n + l)(74n 4 + 148n 3 + 7?i 2 - 67n + 18) (5.54) 

fe=0 

n 

^ #i+ fe - (2n + 1)^+1 " ^ " *Wi + (5.55) 

n , -, 

£ fc^S* = ^ _1) ^ 2) + 2 (2n + l)[fWi - - 2 ( n + !) ( 5 - 56 ) 

fe=0 

k 2 H% = [2H^ +1 - H^} - -(6n 2 +6n+ l)[H 2n+1 - H n ] 

k=o D (5.57) 

+ i(»+l)(3n + l) 
6 

n 1 

]T k^H% = H$- 3 <HM + 3tf (- 2 > [H 2n+1 - H n ] - -n(n + l)(14n + 13) (5.58) 

E** fl S k = fl i- 4) [ 2ff ffl.i- fl i a) ] 

fc=0 

- ~ (30n 4 + 60n 3 + 30n 2 - 1) [# 2n+1 - ff B ] (5 ' 59) 

+ — (n + l)(35n 3 + 47n 2 + 5n - 2) 
60 

n 

Y,k 5 H%l = H^HW + 5H^[H 2n+1 H n ] 

k=o (5.60) 

l)(74n + 67)(3n 2 + 3n- 1) 

oDU 
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H 2n +k = (3n + l)H 3n+1 - 2nH 2n - (n + 1) (5.61) 

fc=0 

™ 1 3 

^2kH 2n+k = --n(3n + l)H 3n+1 +n(2n + l)H 2n + -n(n+l) (5.62) 



fc=0 



™ 1 1 

22 k 2 H 2n+k = -n(2n + l)(3n + l)ff 3 „+i - ^«(2n + l)(4n + l)ff 2n 

fc=o ^ 6 (5.63) 
- -^n(n + l)(40n + 17) 

n 

53 k * H 2n+k = —n 2 {Zn + l)(5n + 3)ff 3n+1 + n 2 (2n + l) 2 H 2n 

fc=o 4 (5.64) 

+ T^n(n + 1)(77ti 2 + 45n - 2) 
48 

n 1 
53 fc 4 tf 2 „+fc = ^n(2n + l)(3n + l)(lln 2 + 5n - l)# 3 n+i 



10 

fc=0 



— n(2n + l)(4n + l)(12n 2 + 6n - 1)# 2 „ ( ' 5,65 ' ) 
15 

1 

1800 



n{n + l)(4692n 3 + 4143n 2 + 467n - 152) 



n 

53 fc 5 if 2 n+fc = --n 2 (3n + l)(14n 3 + 16n 2 + 3n - l)if 3 n+i 



4 

fc=0 



+ in 2 (2n + l) 2 (8n 2 +4n-l)ff 2 „ (5 ' 66) 
+ ^n(n + l)(1036n 4 + 1148n 3 + 177n 2 - 87n + 6) 



Formula (4.7) 



J2 H 2n -k = (2n + l)iWi - nfln - (« + 1) (5-67) 

n 15 
53 fc^2„-fc = n(2n + l)H 2n+l - -n(3n + 1)#„ - -n(n + 1) (5.68) 



fc=0 

™ 1 1 

53 fc 2 iJ 2 „_ fe = -n(2n + l)(4n + l)# 2 „+i - -n(2n + l)(7n + 1)JT„ 

fc=o d (5.69) 

- -U(n+l)(64n + ll) 
3d 

n 

53 k 3 H 2n _ k = n 2 (2n + l) 2 H 2n+1 - -n 2 (3n + l)(5n + 3)Jf n 

fc=o 4 (5.70) 

- —n(n + l)(131n 2 + 51n - 2) 
48 
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" 1 

22 k 4 H 2n _ k = —n(2n + l)(4n + l)(12n 2 + 6n - l)H 2n+1 



k=0 



15 

— n(2n + l)(93n 3 + 66n 2 + 2n - l)H n ( 5 ' 71 ) 
30 

1 n(n + 1) (7932n 3 + 4953n 2 - 43n - 92) 



1800 

" 1 
^ A: 5 i7 2 „_ fc = -n 2 (2n + l) 2 (8n 2 + 4n - 1)^+1 



fe=0 



-n 2 {3n + l)(14n 3 + 16n 2 + 3n - l)# n ^ 5 - 72 ^ 



- + l)(5308n 4 + 4604n 3 + 221n 2 - 251n + 18) 

6 Computer Program 

The Mathematica® [S] program used to compute the expressions is given below. 

Unprotect [Power] ; 
0'0=1; 

Protect [Power] ; 

BernPlus [n_] : =BernoulliB [n] +KroneckerDelta[n, 1] 
PowSum [p_] : = 

Factor [1/ (p+1) Sum [Binomial [p+1 , k] BernPlus [p-k+1] n'k , {k , 1 , p+1}] ] 
HarmFun [p_] : =If [p>0 ,HarmonicNumber [n+1 ,p] -1/ (n+1) "p, PowSum [-p] ] 
HarmSumF [p_ ,m_] : =Simplif y [HarmFun [-p] HarmFun [m] +HarmFun [m-p] 

-1/ (p+1) Sum [Binomial [p+1 , k] BernPlus [p-k+1] HarmFun [m-k] , {k , 1 , p+1}] ] 
HarmSumG [p_ ,m_] : =Simplif y [HarmFun [-p] HarmFun [m] 

+KroneckerDelta [p , 0] HarmFun [m] +1/ (p+1) Sum [ (-1) ~k Binomial [p+1 , k] 

(BernPlus [p-k+1] + (p-k+1) HarmFun [k-p] ) HarmFun [m-k] , {k, 1 ,p+l}] ] 
HarmSumF [p_,m_,s_] : =Simplif y [Sum [(-1) ~k Binomial [p,k] s~k 

(ReplaceAll [HarmSumF [p-k,m] ,n->n+s] 

-ReplaceAll [HarmSumF [p-k,m] ,n->s-l]) ,{k,0,p}]] 
HarmSumG [p_ ,m_ , s_] :=Simplify [ReplaceAll [HarmSumG [p,m] ,n->n+s] 

-Sum [Binomial [p,k] (n+1) " (p-k) ReplaceAll [HarmSumG [k,m] ,n->s-l] , 

{k,0,p}]] 

HarmTable [m_] : =Table [HarmonicNumber [n+1 , i] , {i ,m}] 

HarmTable [m_ , s_Integer] : =Table [HarmonicNumber [n+s+1 , i] ,{i,m}] 

HarmTable [m_,s_] : =Table [HarmonicNumber [s+Mod [i+1 , 2] (n+1) , 

Quotient [i+1, 2]] ,{i.2m}] 
HarmonicSumF [p_ ,m_] : =Module [{t=HarmTable [m] ,u}, 

u=Factor [Coef f icientArrays [HarmSumF [p,m] ,t]] ; 

u[[l]]+Dot[u[[2]] ,t]] 
HarmonicSumG [p_ ,m_] : =Module [{t=HarmTable [m] ,u}, 

u=Factor [Coef f icientArrays [HarmSumG [p,m] ,t]] ; 

u[[l]]+Dot[u[[2]] ,t]] 
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HarmonicSumF[p_,m_,s_] : =Module [{t=HarmTable [m, s] ,u}, 

u=Factor [Coef f icientArrays [HarmSumF [p,m, s] ,t]] ; 

u[[l]]+Dot [u[[2]] ,t]] 
HarmonicSumG[p_,m_,s_] : =Module [{t=HarmTable [m, s] ,u}, 

u=Factor [Coef ficientArrays [HarmSumG [p,m, s] ,t]] ; 

u[[l]]+Dot [u[[2]] ,t]] 

(* Compute some examples *) 
HarmonicSumF [4 , 3] //TraditionalForm 
HarmonicSumF [2 ,2,n] //TraditionalForm 
HarmonicSumF [3,1, 2n] //TraditionalForm 
HarmonicSumG [3 , 2] //TraditionalForm 
HarmonicSumG [3 , 1 ,n] / /TraditionalForm 
HarmonicSumF [2,1,0] //TraditionalForm 
HarmonicSumF [2,1,3] //TraditionalForm 
HarmonicSumF [2 , 1 , s] //TraditionalForm 
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